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In theories with large extra dimensions and TeV-scale gravity, black holes are copiously produced in particle collisions 
at energies well above the Planck scale. I briefly review some recent work on the phenomenology of this process, with 
emphasis on theoretical uncertainties and possible strategies for measuring the number of extra dimensions. 



1. INTRODUCTION 

One of the most surprising and exciting possibilities for new physics at future colliders is the production of 
microscopic black holes Q, Q , which can occur in theories with large extra dimensions 1 Q . String theory suggests 
a scenario in which the Standard Model particles are confined to a submanifold with the familiar three spatial 
dimensions (the 'brane'), while gravity is free to propagate in the full (3 + n)-dimensional space (the 'bulk'). In that 
case the fundamental (4 + Tridimensional Planck scale could be much lower than its effective 4-dimensional value, 
possibly of the order of TeV. Then gravity would look weak on distance scales large compared to the size of the extra 
dimensions, but could become as strong as the other forces at short distances. Parton-parton collisions with centre- 
of-mass energies well above this scale could be treated as (4 + ra)-dimensional classical gravitational interactions. 
Numerical studies suggest that in such circumstances the probability of coalescence to form a microscopic black hole 
would be large jH0- The black hole would be a powerful source of Hawking radiation, decaying rapidly into all 
kinds of Standard Model particles. From a study of the decay products, it might be possible to deduce the number 
and size of the extra dimensions. 

In this paper I shall review some recent work on the above scenario. We shall see that, even if one accepts the basic 
assumption of TeV-scale gravity in extra dimensions, there are many uncertainties and gaps in our understanding, 
which make reliable predictions impossible at present. Nevertheless one can build models and simulations that allow 
one to explore different options for filling in the gaps, and enable experimentalists to start thinking about how they 
might analyse black hole events. 

I apologise for not mentioning all relevant topics and papers in this short review. For more complete discussion 
and references, please see and Chris Harris's thesis 0. 



2. BLACK HOLE PRODUCTION 

On purely dimensional grounds we expect the parton-level cross section for black hole formation to be of the form 

&(s = MIh) = F n nr 2 s . (1) 
where rg is the Schwarzschild radius in (4 + n) dimensions, 

1 



r s = 



J PL 
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and F n is a "formation factor" of order unity. Notice that we assume here that the black hole mass Mbh is equal to the 
full parton centre-of-mass energy this point will be discussed further later. We use the Dimopoulos-Landsberg |3] 



lr rtiis includes warped extra-dimension theories with warp scale large compared to the black hole. 
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definition of the Planck mass, 

M PL = [G(4+n)]^ (3) 

where G(4+„) is the (4 + n)-dimensional Newton constant. For illustrative purposes, we shall usually take Mpl = 1 
TeV in this study. 



2.1. Black Hole Formation Factor 



The formation factor F n in Equation^has been estimated in a variety of ways. A simple geometric argument 
goes as follows. Consider incoming partons that would pass each other at a separation (impact parameter) b = 2r/ l 
where is the horizon radius for a Kerr black hole with the corresponding angular momentum J = b^/s/2 and 
mass -Mbh = s/s. This is assumed to be the maximum impact parameter at which black hole formation could occur. 
Therefore 



a = F n irrl ~ ir(2r h ) 2 
But for a Kerr black hole in (4 + n) dimensions 

r h =r s [l + a 2 ]~^ 

where 

{n + 2)J 



2r h M B n 



and so we obtain 



n + 2 



("geometric"). 



(4) 



(5) 



(6) 



(7) 



This formula, shown by the blue curve in Figure Q follows quite closely the numerical estimate of Yoshino and 
Nambu [3 (green). On the other hand, the latter is only a lower bound on the cross section, obtained by finding a 
closed trapped surface on a particular slice of (4 + n)-dimensional space-time 6]; such a surface must be shielded by 
an event horizon. In this way one obtains a lower bound on the impact parameter for horizon formation, and hence 
on the cross section for black hole formation. 

More recently, Yoshino and Rychkov Q have found a more optimal space-time slice, which leads to a larger lower 
bound on the formation factor, shown in red in Figure 



2.2. Black Hole Cross Section 

Adopting the Yoshino- Rychkov lower bound as an estimate of the black hole formation factor, and assuming again 
that VI = Mbh j we obtain the parton-level cross section, shown as a function of the Planck scale for a 5 TeV black 
hole in Figure [21 For n > 2 extra dimensions, the n-dependence of the formation factor tends to cancel that of the 
Schwarzschild radius, so that the cross section is not strongly dependent on the number of extra dimensions. 2 Values 
of n less than 3 are in any case strongly disfavoured on astrophysical grounds jllT | . 

To estimate the cross section for black hole production at a hadron collider, we must convolve the parton-level cross 
section with the parton distributions in the incident hadrons. The resulting cross section for a pp collider at cm. 
energy 14 TeV (i.e. the LHC) is shown in Figure^ At the LHC design luminosity of 10 34 cm _2 s _1 , this corresponds 
to more than one black hole per second with mass above 5 TeV. 



2 The same is true for the Yoshino-Nambu estimate, only the numerical values are smaller. 
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Figure 1: Yoshino-Rychkov (YR), Yoshino-Nambu (YN) and geometrical (geom) estimates of black hole formation factor. 
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Figure 2: Parton-level black hole cross section. 



2.3. Measuring the Planck Scale 

If we accept the Yoshino-Rychkov result as a reliable estimate of the black hole formation factor, we see from 
Figures |21 and 13 that a measurement of the cross section for a given range of black hole masses would fix the Planck 
mass in a way that is substantially independent of the number of extra dimensions, at least in the astrophysically 
favoured region n > 2. 

Of course, one expects some dramatic changes in the cross section and final state at partonic cm. energies around 
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the Planck scale, due to the onset of strong gravitational scattering. However, to predict those changes one would need 
a quantum theory of gravity, so deducing the Planck mass from them would not be straightforward. Measurements 
well above the Planck scale, on the other hand, can reasonably be interpreted in a classical approximation, as we are 
doing here. 

The problem in any case is to make a reliable measurement of the black hole mass, or more correctly the partonic 
cm. energy for black hole formation. Since we do not observe the colliding partons, this can only be inferred from 
properties of the final state, which will be dominated by the decay of the black hole. 



3. BLACK HOLE DECAY 



Although the formation of a horizon in parton collisions well above the Planck scale seems reliably established, 
the nature and fate of the object thus created is much less clear. The usual working hypothesis has been that the 
evolution of the system has four phases: 

• Balding phase: all 'hair' (characteristics other than mass, charge and angular momentum) and multipole mo- 
ments are lost through gravitational and Hawking radiation, and the object becomes the multidimensional 
generalization of a Kerr-Newman black hole. In fact any residual charge after this phase is probably negligi- 
ble [U, so the Kerr solution is assumed. 

• Spin- down phase: the Kerr black hole loses angular momentum by Hawking radiation and becomes a 
Schwarzschild black hole. 

• Schwarzschild phase: the black hole loses mass through Hawking radiation and its temperature rises until the 
mass and/or temperature reach the Planck scale. 

• Planck phase: the object ('string ball'?) is in the realm of quantum gravity and its fate cannot be predicted. 
It could decay into a few quanta with Planck-scale energies 0, evaporate at the Hagedorn temperature [l^. 
or even form a new kind of stable relic object [l^ . 
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At present, each of these decay phases is subject to great uncertainty. The amount of gravitational radiation emitted 
in the balding phase is of major concern because this constitutes missing energy that would spoil the connection 
between collision energy and black hole mass 3 0, 0, Q ■ In higher dimensions there are solutions to the Einstein 
equation that are not simply generalizations of the four-dimensional solutions, such as black rings Probably 
there are more complicated objects still to be discovered. It is not clear whether such configurations would be 
able to spin down to the Schwarzschild solution, or what their Hawking radiation would look like. Even assuming 
the generalized Kerr solution, the amount and distribution of Hawking radiation during spin-down is still under 
investigation El El El E3 

In the Schwarzschild phase too there are many points that require further clarification. What fraction of the 
Hawking radiation is emitted as detectable Standard Model particles on the brane, and how much escapes into the 
bulk? Is the decay process too rapid for the relationship between black hole mass and temperature to remain valid 
throughout? Do secondary decays significantly distort the Hawking spectrum? And in the final Planck phase, what 
are the consequences of alternative models for the fate of the remnant object? 

Many of these questions await theoretical answers, but some can be illuminated by numerical simulations. For the 
latter, we shall start with the hypotheses that practically all the energy of the parton collision goes into the black 
hole, that the decay is dominated by the Schwarzschild phase, and that the Hawking radiation consists entirely of 
Standard Model particles on the brane. With these assumptions, a detailed picture of the final state can be presented 
and the effects of some of the uncertainties can be investigated. 

3.1. Hawking Spectrum 

With the above assumptions, the spectrum of particles emitted during black hole decay takes the form 

™ oc ^ T"+ 6 (8) 
where as usual the =p applies to bosons and fcrmions, Th is the Hawking temperature 

T H = =±1 oc (9) 

and 7 is a (4 + n)-dimensional grey-body factor H 0, 0, [2^. The latter takes account of the fact that the wave 
function of a particle created in the intense gravitational field near the horizon has to propagate through the curved 
space around the black hole in order for the particle to be observed. By detailed balance, the grey-body factor is 
equal to the absorption coefficient of the black hole for waves incident from infinity. For particles emitted at high 
energies the wavenumber is large compared with the curvature, so grey-body effects are small and the spectrum is 
close to black-body. But at low energy the emission is strongly modified in a way that depends on the spin and the 
number of extra dimensions. 

Figures 0HH taken from Q, show the grey-body factors for scalars, fermions and gauge bosons as functions of the 
particle energy in units of the inverse horizon radius. What is actually plotted here is the absorption cross section 
in units of 7rr|, which tends to 4 at high energies. We see that the main effect in extra dimensions is the suppression 
of low-energy gauge boson emission. 

3.2. Integrated Flux and Lifetime 

Integrating the spectrum in Equation 8, including the grey-body factors, and multiplying by the number of Standard 
Model degrees of freedom for each spin, one obtains Q the total flux of particles of each type in black hole decay. 



s Lower bounds on the black hole mass can be deduced from the trapped surface area QQ. 
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Figure 4: Grey-body factors for scalar emission on the brane from a (4 + n)D black hole. 
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Figure 5: Grey-body factors for fermion emission on the brane from a (4 + n)D black hole. 



Expressed in units of the inverse horizon radius, as shown in Figure the particle fluxes are independent of the 
black hole mass and Planck scale. When the flux in these units exceeds unity, which we see is the case for n > 2, 
the time between emissions is less than the time for a light-signal to travel a distance equal to the horizon radius. In 
these circumstances it is difficult to see how the emission can remain thermal. However, we shall continue to make 
that assumption in the absence of any better understanding. 

Assuming that the Schwarzschild phase of decay is dominant, and that the mass and temperature are relate d by 
Equation El throughout (which we have just seen must be doubtful), the total energy flux can be integrated [ljj 
to find the time at which the entire mass of the black hole has been radiated away. This measure of the lifetime, 
expressed in units of the inverse of the initial mass, is shown in Figure 00 We see that the lifetime falls very steeply 
as a function of the number of dimensions, and indeed can be comparable with the inverse mass when n > 4, even 
for masses well above the Planck scale. When this is the case, the object formed can no longer really be said to have 
an independent existence as a black hole. 
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Figure 6: Grey-body factors for gauge boson emission on the brane from a (4 + n)D black hole. 
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Figure 7: Integrated flux of Hawking emission on the brane in the decay of a (4 + n)D black hole. 

4. EVENT SIMULATION 

4.1. CHARYBDIS Event Generator 

The simulation program CHARYBDIS generates black hole production and decay configurations assuming that 
all the partonic collision energy goes into the mass of the hole, and that the Schwarzschild decay phase dominates. 4 



4 The current public version (1.001) neglects the form factor F n in Equation^ so that the production cross section is somewhat smaller 
than it should be. 
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Figure 8: Mean lifetime of a (4 + n)D black hole. 



Name 


Description 


Values 


Default 


MINMSS 


Minimum mass of black holes (GeV) 


< MAXMSS 


5000.0 


MAXMSS 


Maximum mass of black holes (GeV) 


< cm. energy 


cm energy 


MPLNCK 


Planck mass (GeV) 


< MINMSS 


1000.0 


TOTDIM 


Total number of dimensions (4 + n) 


6-11 


6 


TIMVAR 


Allow Tfj to change with time 


LOGICAL 


.TRUE. 


MSSDEC 


Choice of decay products 


1-3 


3 


GRYBDY 


Include grey-body effects 


LOGICAL 


.TRUE. 


KINCUT 


Use a kinematic cut-off on the decay 


LOGICAL 


.FALSE. 


NBODY 


Number of particles in remnant decay 


2-5 


2 



Table I: Main CHARYBDIS parameters. 




The identities and momenta of the incoming partons and outgoing primary decay products are passed to the HER- 
WIG [24| event generator via the Les Houches interface |^. HERWIG then handles all the QCD parton showering, 
hadronization and secondary decays. 

The main parameters that control the operation of CHARYBDIS are summarized in Table [I] The first four are 
self-evident. The remaining five provide the means to study the effects of some of the uncertainties discussed earlier. 
TIMVAR causes the temperature of the black hole to be updated according to Equation^ after each emission; otherwise 
it is frozen at the initial value. MSSDEC controls whether heavy particles are included in the Hawking radiation: 
MSSDEC=1 allows only light particle emission (up to and including b quarks); MSSDEC=2 includes top quark, W and 
Z emission; MSSDEC=3 includes also Higgs boson emission. 

The parameters KINCUT and NBODY determine how the evolution of the black hole is terminated. If KINCUT= . TRUE . , 
termination occurs when the chosen energy for an emitted particle is ruled out by the kinematics of a two-body decay. 
At this point an isotropic NBODY decay is performed on the black hole remnant. The NBODY particles are chosen 
according to the same probabilities used for the first part of the decay. The selection is then accepted if charge and 
baryon number are conserved, otherwise a new set of particles is picked for the decay. In the alternative termination 
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Figure 9: Primary photon spectrum from black holes with initial masses 5.0-5.5 TeV in n = 2 extra dimensions. Dashed: 
neglecting time-variation of temperature. Dot-dashed: neglecting grey-body factors. Solid: including both. 



Particle type 


Particle emissivity (%) 


GRYBDY= . 


TRUE. 


GRYBDY=. FALSE. 


Generator 


Theory 


Generator 


Theory 


Quarks 


63.9 


61.8 


58.2 


56.5 


Gluons 


11.7 


12.2 


16.9 


16.8 


Charged leptons 


9.4 


10.3 


8.4 


9.4 


Neutrinos 


5.1 


5.2 


4.6 


4.7 


Photon 


1.5 


1.5 


2.1 


2.1 


Z° 


2.6 


2.6 


3.1 


3.1 


W+ and W - 


4.7 


5.3 


5.7 


6.3 


Higgs boson 


1.1 


1.1 


1.0 


1.1 



Table II: Relative numbers of primary emissions from black holes with initial masses 5.0-5.5 TeV in n — 2 extra dimensions.. 



(KINCUT=. FALSE.) particles are emitted according to their Hawking energy spectra until Mbh falls below MPLNCK; 
then an NB0DY decay as described above is performed. Any chosen energies which are kinematically forbidden are 
simply discarded. 

Figure^illustrates the effects of grey-body factors and temperature variation on the energy distribution of primary 
photons. Both tend to harden the spectrum, shifting the peak towards higher energies. They also reduce the total 
number of photons emitted, the grey-body factor by suppressing soft photon emission and the temperature variation 
by reducing the lifetime of the black hole. 

Table [n] shows the relative numbers of primary particles of different types emitted when the parameters TIMVAR, 
GRYBDY, NB0DY and NB0DY have their default values, compared to the values obtained by integrating the theoretical 
spectra. The minor discrepancies are due mainly to kinematic constraints, particle masses and charge conservation. 
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Figure 10: Missing transverse energy for various processes at the LHC. 

4.2. Event Characteristics 

Turning from single-particle spectra to overall event characteristics, one interesting feature is that black hole 
decay is associated with large missing transverse energy, due to copious primary and secondary neutrino emission. 
Figure ITU1 fro m [al , shows a comparison with expected QCD and SUSY missing E T distributions generated using 
HERWIG |3i EU- We see that the black hole missing Et is typically larger even than that in supersymmetric 
processes, where missing energy is mostly carried off by a pair of neutralinos. The effect is partly due to the larger 
mass of the black hole, relative to the assumed SUSY scale. However, a large cross section with large missing energy 
would clearly be a good initial indicator of black hole production. 

4.3. Measuring Black Hole Masses 

The large missing energy in black hole decay poses a problem for the reconstruction of the mass of the black hole 
from its decay products. In Q we found that a cut on missing Et < 100 GeV was necessary for a useful mass 
resolution of around 4%, i.e. ±200 GeV at 5 TeV, as illustrated in Figure HP At this low value of missing Et, QCD 
background has to be controlled by requiring at least 4 high-i/r jets. 

5. DETERMINING THE NUMBER OF EXTRA DIMENSIONS 
5.1. Fitting Emission Spectra 

In the event that black holes are produced at the LHC, the quantity of principal interest will be the number of 
extra dimensions, n. Given the sensitivity of the Hawking temperature to n in Equation^ it would seem that fitting 
the emission spectra for various black hole masses would allow one to extract this quantity 0. However, we found 
in that the theoretical uncertainties, together with the distortion of the spectrum by secondary decays, would 
make it difficult to have confidence in such a measurement. 
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Figure 11: Black hole mass resolution after cuts. 
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Figure 12: Fits to Tjj vs A/bh from primary electron spectrum for n = 2. (a) TTMVAR off; (b) TIMVAR on. 



For example, Figure El shows the effect of possible time variation of the Hawking temperature on a fit to the 
primary electron spectrum (assuming this could be unfolded cleanly from the data). For events generated with the 
temperature frozen (left), the fit gives a result consistent with the input value n = 2, whereas re-thermalization 
between every emission (right) systematically shifts the fit to higher values, due to the higher average temperature. 
Since the true situation would presumably lie between these extremes, the true value of n could not be extracted 
without a deeper understanding of the decay process. 

Figure ^5] illustrates another difficulty in fitting the emission spectrum as a function of black hole mass, this 
time due to kinematic effects at higher energies. The limitation to emission energies less than half the total mass 
significantly truncates the spectrum generated according to Equation |H1 at low masses and/or large values of n. 



5.2. A Possible Observable 

In order to avoid the low-energy region of Hawking emission, where secondary decays distort the spectrum, and 
the highest energies where the kinematic cutoff takes effect, we examined the region of high but not extreme 
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(a) n = 2. (b) n = 4. 

Figure 13: Energy distribution of primary decay products vs Mbhi showing the effect of the kinematic limit E < AfBH/2- 
The colour code indicates the number of particles in each bin per 1000 events. 



energies, E ~ _E cut — Mbh/2 — Ed where Ed is a few hundred GeV. Particles in this region are also less sensitive 
to time variation of the temperature, since they tend to be emitted early in the evolution of the black hole. This 
tendency is further enhanced by demanding that the highest-energy emission should have energy E m&K > E cut . We 
therefore looked at the fraction F of events satisfing this cut as a function of Mbh- 

As shown in Figure for Ed = 400 GeV this observable is indeed relatively insensitive to the uncertainties 
represented by the CHARYBD1S parameters TIMVAR, KINCUT and NB0DY. The upper and lower bounds were obtained 
by integrating the Planck spectrum from E cn t up to infinity and up to Mbh/2, respectively, and include a ±200 
GeV systematic uncertainty in the black hole mass. The value of F always lies within, or at least is consistent with, 
the expected band. Therefore a fit to this observable should have less model dependence than analyses based, for 
example, on a fit to the full emission spectrum. 

In Q we combined a fit to this observable and a cross-section measurement with an assumed error of 20%, to make 
a joint determination of the Planck mass and the number of extra dimensions. As shown in Figure ITBI for Mpl = 1 
TeV and n — 4 the method gives an unbiased estimate of these quantities, with la errors of AMpl/A/pl ~ 15%, 
An ~ 0.75 (strongly correlated). 



6. CONCLUSIONS 



If there are indeed extra dimensions large and numerous enough to reduce the fundamental Planck scale to the 
TeV range, black hole production at accelerators, with substantial cross sections, is a real possibility. Such a 
discovery would arguably be the most profound advance in fundamental physics since general relativity. It would 
also constitute the first observation of Hawking radiation. However, as I have repeatedly emphasised, there remain 
many uncertainties about the precise nature, evolution and ultimate fate of the objects that would be formed. The 
strategy adopted here has been to start with a crude but flexible simulation framework, and to look for observables 
that are less sensitive to some of these uncertainties, in order to stimulate thinking about how the fundamental 
parameters could be measured. If real data on this process do become available, there will undoubtedly be an 
explosion of theoretical and experimental activity that will rapidly reduce the uncertainties. 
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(b) No time variation. 




5000 



6000 



7000 



8000 



9000 10000 
M B „ (GeV) 




5000 



6000 



7000 



8000 



9000 



10000 

M B „ (GeV) 



(c) Kinematic cut on. 
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Figure 14: Fraction of events with E max > Mbh/2 — 400 GeV vs Mbh for different options, all with n = 4. Upper and lower 
bounds explained in the text are also shown. 
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Figure 15: Combined measurement of Planck mass and number of extra dimensions. 
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